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Abstract 

We study the Godel universe through worldlines associated with mo- 
tion at constant speed and constant acceleration orthogonal to the instan- 
taneous velocity (WSAs). We show that these worldlines can be used to 
access every region — both spatial and temporal — of the space-time. 
We capture the insights they accord in a series of sketches, which extend 
significantly the Hawking and Ellis picture of the Godel universe. 

1 Introduction 

The Godel Universe jT] has several remarkable properties. It is homogeneous 
(manifestly so in all studied versions of 'rectangular' coordinates) rotationally 
symmetric (quite obviously, when viewed from an axisymmetric cylindrical sys- 
tem), anti-Machian (owing to its rotation as a rigid body in an inertial frame) 
and acausal (because its light-cone structure gives rise to closed null and closed 
time-like curves (CNCs/CTCs)). 

Its geodesies [21 E] inevitably reflect one or more of its underlying properties. 
In both the cylindrical and Kundt's version of rectangular coordinates, they 
project onto circles in the (r, (j)) and (x, y) planes respectively (their closure being 
a direct consequence of the universal rotation) and correspond, quite generally, 
to motion at a constant speed, v. As v increases, the radius of the circle increases 
as well, but remarkably, the time taken to circumnavigate it completely, as 
measured by an observer situated at any one of its points, decreases! The largest 
geodesies are null curves, stretching out to the Godel horizon which marks the 
boundary of causal behaviour as also of the region which the origin can access 
through geodesies. 

By subjecting our projectiles to a uniform acceleration. A, orthogonal to 
their instantaneous velocities, we can force them onto orbits of ever larger radii 
and can thereby bring them back to the base-station ever earlier. So much so 
that for a critical acceleration, Actc{''j)j we could bring them back instantly, 
and for an even larger A, could bring them back before they were launched! 

All of these properties can be depicted diagrammatically and some have in 
fact already been so captured. In an enlightening picture of the Godel universe 
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in cylindrical coordinates, Hawking and Ellis [4 brought out beautifully the re- 
focussing of photons diverging uniformly from an arbitrarily chosen origin, the 
symmetrical tilting of light-cones in the direction of rotation, the emergence of 
CNCs at the Godel horizon and their turning into CTCs thereafter. Pfarr[S] 
took the first step towards developing a similar sketch in Kundt's coordinates, 
by tracing the corresponding geodesies, CNCs and CTCs along with their planar 
projections. 

In this paper, we take Pfarr's discussion to its logical conclusion. We start by 
showing that the extensions of geodesies he considered, correspond to motion at 
uniform speed and magnitude of acceleration, orthogonal to the instantaneous 
velocity. The WSAs, as we shall refer to these extensions, have all the symme- 
tries of geodesies and in particular the property of refocussing, which now gives 
rise to a truly intriguing series of diagrams: Indeed, as we approach Actc{v), 
the points of divergence and re-convergence draw ever closer, and finally coin- 
cide as A = AcTc{v)- As we increase A further, we eventually end up with 
open trajectories, and for even higher A, trajectories which curve oppositely to 
the geodesies. Trajectories for these values of A have not, as yet, been studied 
in the literature. In fact, WSAs provide a very convenient and insightful way of 
studying all the interesting phenomena associated with the Godel universe. 

We also present ways of visualising the general mechanism for the emergence 
of CNCs/CTCs. Indeed, by explicitly drawing a set of representative light- 
cones along the latter, we show how the tilting of these cones along with their 
opening out to include a spatial axis (in addition to the temporal one) allows 
us to complete the CNC/CTC. We first show this for cylindrical coordinates 
and then, in the backdrop of a manifest translational symmetry, for Kundt's 
coordinates. 

Finally, we note that v and A are the control variables so to say for navi- 
gating Godelian space and time. By supplying the formulae and a diagram to 
show how these can be adjusted to access a coordinate time of choice, within 
a specified interval of proper time, we arrive at a practical travel guide for the 
Godel Universe. 

The plan of the paper is as follows: In section [21 we set up the equations of 
motion for planar geodesies passing through the origin in cylindrical coordinates 
and show that the geodesies have a constant speed. Further, the geodesic of 
speed V projects onto a circle in the R-(j) plane of diameter D = v/^/2. In section 
12. H we obtain a family of worldlines from the geodesies by stipulating that the 
diameter, D, of the trajectory be independent of its speed, v. We argue that the 
worldlines so obtained correspond to motion at constant speed and magnitude 
of acceleration. We calculate this acceleration by comparing each accelerated 
worldline with a geodesic which shares the same trajectory in the R-cf) plane. 
In section [2.2) we determine the coordinate time, t, and proper time, s, taken 
to reach a given point on the trajectory of a worldline. We use these results 
in section 12.31 to list several noteworthy properties of these worldlines, which 
we additionally capture diagramatically in section 12.41 One of these diagrams 
lays bare an intuitive mechanism for the emergence of CTCs, which generalizes, 
in a certain sense, the one captured by Hawking and Ellis. In section [31 we 
rederive several of the above results in Kundt's coordinates. We moreover use 
the manifest homogenity of this coordinate system to remove from the orbits 
of section [21 the constraint that they must pass through the origin. In section 
13.11 we mimic the results of section 1^21 in Kundt's coordinates. In section [221 
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we visualize these results from several perspectives. Finally, in section 31 we 
summarize our findings. 



2 Kinematics - Cylindrical Coordinates 

Our first objective is to construct 'planar' wordlines in the Godel universe (i.e. 
worldlines on the plane orthogonal to the z-axis) which start from the origin 
and travel with a constant speed and magnitude of acceleration. Such worldlines 
are considered because they fully exploit the symmetry of the universe, they are 
reasonably simple to handle and they are sufficient to illustrate the strange 
properties of the Godel universe. 

We first construct these worldlines in the cylindrical coordinate system (t, r, (jj) 
[T], for which the line element is: 

ds^ = (dt + 2 sinh^ f—\ - ( dr'^ + sinh^ r cosh^ f (1) 



where f — ct;r/-\/2 and uj is the angular velocity of rotation of the universe. 
The dz'^ term has been dropped as we consider only trajectories restricted to 
the plane of rotation of the universe (the r-(j) plane) . The calculation simplifies 
considerably |6j, if we rewrite the metric in terms of i? = tanh(f) (which maps 
the entire space-time onto the interior of the unit circle): 

The effective Lagrangian for a particle moving geodesically is then 

?2 I n2l2 



C - g^^x^x- -\^+ ^(^_^2) j - ^2(1 _ ^2)2 - ' (3) 

where e = 0/1 for photons/material particles and i — dx/dX. For material 
particles, X — s. 

Since C is independent of t and (j), there are two constants of motion: 

(4) 




(5) 

Note that = for geodesies passing through the origin, while pt is related 
to the speed of the geodesic. Indeed, the speed of a moving particle as seen by 
comoving observers is given by [7] 

2 {goidx'Y - googtjdx'dx' g^ui'^i" 
= 1 nr^ Tlyi = 1-5007 ^5- (6) 

where j, j index only the space coordinates. In the present case, using 500 = li 
^ and ([4]), we get a constant 



t + 



= 1 - — T = 1 - -2 (7) 
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Figure 1: The i?(/)-projections of the geodesies for w = 1 and v = l/\/2. The 
Godel Horizon corresponds to R— l/-\/2 and spatial infinity to i? = l.Note that 
in these coordinates, this plane is isotropic: particles diverging from the origin 
with the same speed describe circles of the same radii. 

Eliminating t between ^ and ([5]) to obtain (j) = io{l — E?)pt and substituting 
this and (g]) into ^ to obtain 7? = cj(l - R'^)pt^^{l - e/p1) - R? , we get 

which readily integrates to 

i? = sin((/. - (/.o) (9) 

It follows that geodesies which pass through the origin are circles in the 
R-(j) plane with diameter D = u/\/2- Solution ^ corresponds to a projectile 
from the origin with initial speed v in the (pQ direction. The maximum radial 
coordinate it attains is Rmax = v/\/2, ie. rmax = iV^/'-^) tanh^^(u/-\/2), which 
for photons becomes r^^^ = (\/2/w) tanh"^(l/V^) = (\/2/w) sinh"^ 1. No 
geodesic from the origin can go beyond the circle r = r^^x^ the Godel Horizon. 

2.1 Worldlines of Constant v and A 

To see how geodesies can be extended to worldlines with a constant acceleration, 
it is instructive to examine them in the non-relativistic limit. For w <C 1, the 
geodesic trajectories given by © are restricted to i? ^ 1. In this hmit, R ~ 
tanh~^ R = rLu/\/2 and hence the trajectories in the r-(j) plane are approximately 
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r = (v/lu) sin(0 — (po), that is, circles of diameter d = v/lu. This result can be 
obtained from a different perspective: 

The comoving frame of reference (dt, dr, d^) is not inertial. At any given 
point, it looks like a rotating frame of angular velocity to, with the axis of ro- 
tation passing through that point There is a Coriolis-like inertial acceleration 
2ujv on a particle moving with non-relativistic speed v, perpendicular to the di- 
rection of motion [5]. This provides the centripetal acceleration 2v'^/d required 
to keep the particle moving with speed v on the circle of diameter d. Setting 
2lov = 2v'^/d and we get d = v/uj. Of course, as seen from the local inertial 
frames, the geodesic has zero acceleration and is locally 'straight'. 

Now consider the motion of a particle with constant speed v <^ 1 along a 
circle of diameter d ^ v/uj. Clearly the required centripetal acceleration 2v'^/d 
no longer equals the Coriolis acceleration 2ujv. An additional acceleration A = 
2ojv — 2v^/d in the radially outward direction would be required to maintain this 
trajectory at this speed. As seen from the local inertial frames, these worldlines 
would be accelerated and would have just this acceleration A ~ 2ljv{1 — v/ujd). 
These worldlines have all the symmetry of the geodesies they have been derived 
from. All points on them are equivalent. They have a constant speed and an 
acceleration of constant magnitude, perpendicular to the velocity. All we need 
to do is extend this construction to the relativistic case. 

This can be done by considerng the worldine of a particle W moving at 
constant speed v along the trajectory 

R = D sm{(l) - (j)o) (10) 

in the R-(j) plane, with D specified independently of v. The worldline is a 
geodesic only if w = Dy/2. In general, each worldline of this 2-parameter family 
(labelled by v, D) has an acceleration, A{v, D), which can be calculated neatly 
by comparing its worldline with that of a geodesic G of speed Vg = D\/2 which 
projects to the same trajectory in the R-(f) plane. 
The 4-acceleration of the W worldline, a;''(s), is 



ds^ ds ds 



while that of the geodesic G, 2^g(s), is 

d'^xP dxP dx'^ 

0=^+rLT^T^ (12) 
ds^g ^ dsg dsg 

The worldlines of W and G differ only in the t and s at which they reach 
a given point on the R-cf) trajectory. So if we move along the two worldlines 
maintaining correspondence between points of same {R,(j)), then dRg = dR, 
d(j)g — d(j), the metric and FJ^j, (which depend only on R) are the same at 
the corresponding points of W and G, but dsg and dtg differ from ds and dt. 
(Subscript g labels quantities for the geodesic G). 



^ Every point can appear to be on the axis of rotation only in a curved spacetime. Con- 
ventionally, the non-relativistic limit of the Godel solution is taken by fixing a rotation axis, 
and introducing a gravitational field to exactly counter the centrifugal acceleration[8]. In any 
case, the inertial acceleration at a given point has only the Coriolis term. 
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The latter can be related by rearranging the defintions of v and Vg 



2_ 2(di?2 +i?2^(/,2) 

~ ((1 - R'^)ujdt + 2R^d(l)f 

2 _ 2{dR^ + R^dcj)^) 
^ ((1 - R^)Ludtg + 2R^d(l)f 
to the form 

ds 

dSg 

dt 



2{dR^ + R^d^^) 



((1 



R^)udsy ^ 

2{dR^ 



2(di?2 H 

- R^dcjP ) 



R^d(P^) 



((1 - R^)udsgy + 2(di?2 + i?2d02) 




_ dig ^ /i^ _ \ 1 

C?Sg dSg \ V ) 



(13) 
(14) 

(15) 
(16) 



Since {ds/dsg) is constant, d'^x^ /ds^ — {ds g / ds)"^ {d"^ x'' / ds"^) and since {dt/dsg)~ 
(dtg/dsg) is likewise a constant, while R and are identical at corresponding 
points, d'^xP [ds^g = d'^x^/ds'^g. Therefore using ITTI) reduces to 



ds 



ds2 



dxf" dx" 

dSg dSg 



= 1^1 re, 

ds 



da;'' dx" 

dSg dSg 



dx^ dx^g 

dSg dSg 



(17) 



Next, since dx^/dsg differs from dx'i/dsg only in the ^-component and 



Fqq = for a stationary metric with constant 500, only Fq^ (^^^f ^'JT' 

contributes, and in view of (fTOjl gives 

2 




(18) 



Substituting into this, the only non-zero Fq^ viz. T^^ = 2R/ (1 — R^) 
-oj/R, we getd 



and Fq;^ 



= ujR 



V 



2uj ( ds 



ds 



^_^2 Vw(l-i?2)ds 



2R dR ^d(j) -ldR\ 

'^7^' (19) 



dso 



R dSg J 



Using (fTS]) to substitute for (dsg/ds) and (fT4| to replace the dsg in (dR/dsg) 
and (d(f)/dsg) in terms of -y/ di?2 R^dcf)'^ , simplifying the obtained fractions in 
dR and d(/i using dR ~ dcjy^/ D"^ — i?2 and substituting = D\/2, wc find that 
the contravariant components of the 4-acceleration are 



AP = A 



Lojl^R^) 
V2 



( 2i?x/D2 - j?2 ^ ^/p-^-R-^ ^ 



(20) 



its covariant components are 



A„ = A 



V2 



Cj(l-i?2) 



0, 



^ R 
15' D 



\/D2 _ J12 



(21) 



^The components are listed in the order {t,R, 
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and its magnitude is 



Thus the worldhnes have a constant magnitude of acceleration A determined 
entirely by the constant speed v and the diameter, Z?, of the trajectory. In 
general u^A^j^ = 0, and in this case Aq = Q. Thus uM^ = 0, which implies that 
the 3-acceleration is perpendicular to the 3-velocity, as desired. 

Note that A is signed. When ^ > 0, the acceleration is opposite to the 
inertial Coriolis acceleration!^ This can be verified most easily by taking the 
t; ^ 1 hmit. We get A « 2ivv{l—v/ujd), where d is the diameter of the trajectory 
in the r-cf) plane (i.e. Lud/y/2 := tanh~^ D « D). This result agrees perfectly 
with the non-relativistic arguments given earlier. 

The magnitude of the 4-acceleration, A, is equal to the magnitude of the 
3-acceleration in a Local Inertial Frame (LIF) instantaneously at rest with W. 
Indeed, the thrust that a rocket of rest mass mo must generate to travel along 
the worldline W (as seen in its own rest frame) is itiqA, directed perpendicular 
to its velocity with respect to the comoving frame. In particular, if this thrust 
is generated by ejecting gas/accelerated ions at relative speed Vreh the rate of 
ejection of mass (in the rocket's rest frame) is given by 



drriQ \ 



/I _ ,,2 

^ '■''^moA (23) 



ds J Vrel 

The LIF is however not meaningful for a photon and indeed A diverges as 
V —> 1. In this case, the required acceleration may be generated through 
forces applied by appropriately stationed objects comoving with the universe. 
The 3-force to be applied at any given point can be determined by transforming 
to the comoving LIF the 3-force, i^est — wo^: we have been dealing with in 
the rest LIF of W . This can, in turn, be done by applying to ttiqA a Lorentz 
boost of speed v perpendicular to iVe^t . The result is 



FroLv^na = moAVT^=EA{l-v' 



V 



2™i5|l-^) ,24) 



where E = mo / Vl — v'^ is the total energy of the particle that is to be forced 
along a WSA of speed v and diameter D. Note that F^omoving perfectly well- 
behaved and hence meaningful at w = 1. As a concrete example, we could have 
an electromagnetic pulse {v = 1) of energy E relayed across a series {N ^ 1 
of comoving stations positioned along a circle of diameter D on the R-cj) plane 
^comoving then approximately the impulse imparted to the photons at each 
station, divided by the transit time between successive stations. 



■^This means that the acceleration is 'outwards' to the trajectory, but only as long as 
the acceleration is smaller than the Coriolis term. When the acceleration A > exceeds 
the Coriolis term, the sense of the trajectory reverses and the acceleration becomes inwards. 
These trajectories can be described by _D < and are discussed in Section 12.31 

^This is only to be expected because A is the magnitude of iPxi^ /ds^ , and ds = for v = 1. 

^With D > DcTc\v—l = 2v^/3, this could be used for sending signals to the past. 
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Before analysing these results further, we derive the expressions for the co- 
ordinate time t, and the proper time s taken to get to a point on the worldline's 
trajectory in the R-(j) plane. 



2.2 The t and s integrals 

We start by rearranging (|13p . after eliminating R from it using R ~ Dsm{(f> — 
4>q), to the form 

dt _ DV2/v - sin^icj} - (j)o) 
^ 1 -D2 sin^ict)- (j)o) 

ds DV2 

ui — = ^ (25) 

d<j) V 1- D^sm^{(j)-^o) 

These can be recast to 

dtanid} — (hn) 



uj{t - to) 



2d(j)- {2- DV2/v) 



1 + (1 - D^)ta.ii'^{(f)- 



DV2 r /-^ dtan(0-0o) 

V 1 + (1 - £12) tan^((/« - (/)o) 

and readily integrated for \D\ < 1 to[^ 

uj{t — to) — 2(0 — 0o) — ^ dV2/v g^^^^g^^ i\/\ — D'^ tan((/) — (/jp)') 

So) = ~ - J^^^^—j arctan [\/\ - D"^ tan(0 ~ 0o)^ (27) 



for ID] = 1 to 

^{t-to) = 2(0- (/«o) - (2 - \/2/v) tan( 



uj{s~so) = V2tan((/)- 0o) (28) 

V 

and for \D\ > 1 to 

uj{t~to) = 2(0 - 0o) - ^ ~ £'\/2/ t; ^^^^^-1 /'^j)2 „ I tan(0 - ^p)") 

cj(s-so) = ^ VIp^ ^^"^'^ (^^' ~ ^ ^^''^'^ " '^"^) ^^^^ 

Now that we have t, R and s as functions of 0, we have a complete parametric 
representation of the worldlines. All the above results have been derived with a 
worldline with 13 > in mind. But the results are valid for D < as well, if we 
consider negative D to correspond to trajectories of diameter \D\, but traversed 
in the opposite sense. Q 



^Here the arctan has been defined with range [0,7r), and arctan(fc tan(7r/2)) := tt/2. This 
makes this solution vaUd upto (/) = 4>0 + ""i ie- up until the worldline has completed one orbit. 
The solution can be extended to any number of orbits, by adding mr to the arctan, where n 
is the number of completed orbits. 

^For the D < case however, ip decreases with increasing proper time, and the worldline 
that leaves the origin with </< = 0o returns to it with = </<q — tt. Hence the arctan used in 
I I27I I for D < must be defined with range (— tt, 0], along with arctan{/c tan(— 7r/2)) := — 7r/2 
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Figure 2: WSA's of a given speed v. This diagram summarizes the relation- 
ship between Z?, A and the projection of the trajectories of the correponding 
worldhnes in the R-cf) plane, for a given speed V > l/y/2. (If v < l/V^, the 
diagram would be the same except that there would be no CTC, as Dqtc > 1 
for V < I/V2.) 

2.3 The Properties of the WSAs 

As we have seen, WSAs depend on two parameters, v E [0, 1] and D e {—00, 00). 
D is the diameter of the WSA's i?-0-projection. It is positive for orbits coro- 
tating with geodesies, and negative for those counter-rotating. The broad be- 
haviour of WSAs for a fixed v and variable D is depicted in Figure [21 

Note that the geodesic of speed v has a fixed diameter D = in the 

R-(p plane and, of course, ^ = 0. For an inward acceleration ^ < 0, the applied 
force assists the inertial Coriolis acceleration and the trajectory shrinks to a 
D < v/^/2. As the inward acceleration A — > ~oo, D ^ and the trajectory 
shrinks to a point. For an outward acceleration A > 0, the applied force opposes 
the Coriolis term, and makes D > v/^/2. The trajectory continues to expand 
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with increasing A, upto A — > Aopen+, the critical 'Escape Acceleration' 

At A — Aopen+, D ^ I, and we get the critical escape trajectory, namely, a 
semicircle with the point at i? = 1 removed. As we increase the acceleration 
further, D increases, and the trajectories turn into minor arcs starting from 
R = and terminating at i? = 1. The proper time taken to approach the R = 1 
point is infinite for all of them, as we would expect. 
When the acceleration reaches the value 

2VU! , . 

Astraight = ~^ ^ 

it exactly cancels the Coriolis term, and the trajectory becomes a straight radial 
line from R = to R = 1. This corresponds to I? ^ cx3 (as also to D ^ — oo). 
Higher accelerations produce circular arcs curving in the opposite direction, 
and hence corresponding to D < 0. Note that the acceleration is now inwards, 
while the Coriolis force is outwards. As A increases beyond Astraight, D starts 
increasing from — oo, and the trajectories are circular arcs of the opposite sense 
terminating at i? = 1. This continues till A reaches the value 

Aope„- = ^(l + ^) (32) 

for which D ~ —\. This corresponds to the critical escape trajectory in the 
opposite sense — likewise a semicircle terminating at i? = 1. Thus for A G 
[Aopen+, Aopen-], wc get \D\ > 1 and hence, open (unbounded) trajectories. 

For accelerations higher than Aopen-, the trajectories are circles in the R-(f) 
plane of diameters D E (—1,0), directed oppositely to the spacetime's geodesies. 

The proper and coordinate times required to circumnavigate the closed orbits 
identified above have a series of interesting properties. For D < 1, the worldline 
which leaves the origin along (j) = (pQ, returns to it along (j) ~ (f>o + n. Setting 
4> = 4>o + TT in (|27p we obtain times At and As required to complete the orbit, 
as measured by the observer O stationed at the origin, and the observer W 
traveling along the worldline, respectively. 



DV2/v^ 



'^^r'ound-trip ~ , , ( ^ — I ('^'^) 




(34) 

For geodesies, D = i?/\/2 and the above expressions reduce to 

A .geodesic 

round— trip 



/. geodesic _ n \/ ± v 

'-^"round-trip ~ . . r. TTTT- V'^'^'' 
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For w — > 0, both O and W observe the time for the round trip to be tt/w. 
Geodesic round-trips of higher v complete sooner, according to both the ob- 
servers. For nuU geodesies, that is at w = 1, the time observed by W is as 
expected, while that observed by O is (7r/ijj)(2 — \/2). 

For a fixed v, Asround-trip is a monotonically increasing function of D, which 
means that spatially larger orbits take longer to complete, as observed by the 
traveler W. At small D, Asround-trip varies linearly with D and as 13 — > 1, 
Asround-trip oo which simply means that orbits of infinite spatial extent 
take infinite proper time to complete. Hence the behaviour of Asround-trip is 
perfectly normal. 

By contrast, the behaviour of Atround-trip is completely normal only at 
small D where it varies linearly with D. Indeed, Atround-trip is not monotonic 
in D for fixed v. It has a maximum at 13 = 1 /v^/2. As long as w < 1/ v^, it does 
not behave strangely for any D £ (0, 1) and as 13 — > 1, Atround-trip +oo. 
But for V > l/v^, as Z3 — > 1, Atround-trip —>■ — oo! This means that by moving 
along a sufficiently large orbit at a speed greater than W can travel as 

far into the past (as measured by O) as desired. In particular, W can choose to 
move along a Closed Time-like Curve for which Atround-trip — and hence 

For any v > l/\/2, worldlines with D > Dctc or equivalently A > Acre are 
Past- Travelling Worldlines. Of course D must be < 1 for the orbit to be closed 
in the first placed 

As noted above, there is a second set of closed orbits, directed oppositely to 
the space-time's geodesies. For this set, the higher the acceleration, the smaller 
the diameter, to the point that D 0— a.s A ^ +oo. We can calculate the 
time taken to complete any of these orbits exactly as in section 12.21 The only 
changes are that D < and our integral over goes from 0o to 0o ~ ^r, since 
now decreases with increasing proper time. The results are: 



DV2/1 



^^round-trip ~ , . I P. — ^ I (^''') 



= VT^^^ (38) 



Vl-i32 



In fact, these D < formulae can be obtained from their Z3 > counterparts, 
((55)1 and ([331), simply by replacing s ^ —s and t ^ —t. Further, for these orbits 

rAt^<o = ^ (As^<'> \+—( ^ - 1 

\^ round-trip) V^l~^ \ round-trip) ^ \Vl- D'^ 

Hence (At^olild-tnp) > {^s^ould-tnp) > 0- In particular, there are 

no CTCs of the reverse sense. Moreover, roughly speaking, the curvature effects 

^From the viewpoint of fuel efficiency, the optimum CTC corresponds to Vopt = 



•\/ V3/2 0.9306. Indeed, by (I23I I. the mass ejected over the round-trip is Am = 
jtIq (1 — exp (^—AAsround-trip^/^ ~ ''Jig; /^rei) ) • Clearly, the minimum of Am coincides 
with that AAsround—trip- ^opt minimises the latter expressed as a function of v alone using 
| [22t . Il34t and 
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aid the time-dilation effect i.e. W can 'travel further into the future' of O than 
possible in a flat spacetime, by making a round trip at a given speed v. 

All these results are best summarized by a Atround-trip-^Sround-trip dia- 
gram (See Figure[3|), which we draw for convenience in terms of the dimensionless 
round-trip times - T := LoAUound-trtph and S := LoAsround-trip/T^- 

There are 3 distinct regions visible on this T-S plane, differing in the number 
and nature of solutions of (w, D) for each point (T, S) in them. □ 

The first region, S € (0, oo) and T G (—00, 5), is bounded by the hue T — S, 
which corresponds to the non-relativistic limit, d — > and Z? — s- 0. In this region, 
there is exactly one value of v and D, with D > 0, for each value of T and 5* 
namely: 



D=J1-1„ v=(l 



F2' V 2(^2-1) 



-1/2 



where 



F = F{S, T) = (2 - T) + -i/ i [(2 - T)2 + - 2] 



2 

Note that W can travel as far back into the past of O as desired, with as small a 
proper time lapse as desired, by making an appropriate choice of v and D (that 
is, bringing them both sufficiently close to 1). 

The second region corresponds to 5 € (0, 1) and T G {S,2 — \J2 — 5*2), lying 
between the above mentioned line and a circle of radius \f2 centred at T = 2, 
S' = 0. In this region, each value of T and S gives rise to exactly three solutions 
for V and D. Two of the solutions correspond to Z? > 0, and are given exactly 
as in the first region, except that in this case there are two possibilities for F: 



F^F±iS,T)^{2-T)±^^ [(2 - T)2 + 52-2] 
The third solution corresponds to < 0, and is given by: 



G2' V 2(G2 - 1) 

where 



G = GiS, T) = (2 + T) - y i [(2 + T)2 + S^-2] 

The rest of the T-S plane constitutes the third region, where there is exactly 
one solution for v and D, with D < 0. This solution is exactly the same as the 
D < solution in the second region. 

These results can be easily visualised by imagining the iso-w and iso-D lines 
on the diagram as two sheets, a sheet for Z? < and a 'folded' sheet for D > 0. 
There are as many solutions of {v, D) at any (T, S) point as the total number 
of overlapping surfaces of both these sheets at that point. 

We end this section by pointing out an interesting property of low velocity 
geodesies. First note that 



^^round—trip "^/l"" 

^ ^^round—trip 



^The solutions {v,D) for each {T,S) can be obtained by inverting 1331 1. 1134 l l and 11371 . 
I|38|l. The equations are quadratics, and the solutions must be checked by back-substituting, 
to discard spurious solutions. 
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Figure 3: Lines of constant v and constant _D on a T vs S* plot. Each set of lines 
can be visualised as a grid drawn on a stretched sheet. There are two distinct 
sheets, one for D e (—1,0) which stretches over the entire region T > 5 > 
(drawn in light gray). The other sheet for D € (0, 1) (drawn in black) has a clear 
fold along the arc T = 2 — ^2 — S'^, S £ (0, 1). A single fold of this sheet covers 
the entire T < S, S > region, while two folds of this sheet cover the region 
between the line T = 5 and the arc T = 2 - \/2 - S^, S G (0, 1). First look 
at the D > sheet. The lines of constant i> for i> < l/-\/2 always have positive 
slope, and T always increases with increasing S. However for v > l/\/2, there 
is a maximum beyond which T begins to decrease with increasing S. There is a 
finite S at which T = 0, and for S +oo, T — oo. The lines of constant D 
always have positive slope - at fixed D if traveler W observes a longer time, so 
will the base station O. For D > l/\/2. the lines of higher D lie entirely below 
lines of smaller D. So once W crosses the Godel Horizon, the further W goes 
keeping S fixed, the smaller is T. In the D < sheet, both T and S increase 
monotonically with increasing v as well as with increasing \D\. 
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R=l 



Figure 4: A map of the future lightcones of the Godel universe in the cylindrical 
{t, R, (j>) < R < 1 coordinates. All the light-cones contain the future-directed 
<-axis. Note that the Godel horizon is a CNC. 

This is just the usual time dilation effect in a twin-paradox-like scenario, as the 
effect of the spacetime's curvature vanishes in the limit D 0. The 'traveller' W 
measures a time which is smaller by a factor of \/l — v'^ than the time measured 
by the 'stationary' observer O. Note that in this case O is moving geodesically, 
while W is accelerated. 

Contrast this with a scenario where W is also moving geodesically, but with 
a non-relativistic speed w ^ 1. Using (f35|) . we see that 

^^qeodesic \/\-~V^ . a. 

-—^ = — , = 1 + Oiv ) 

^tgeodestc 2^1-^2/2-1 

Since W makes a round trip in this case and returns to the base station O, 
without either of them ever accelerating, they are both on an equal footing. 
Indeed they measure the same time for the journey upto 0{v'^). At this order, 
the effect of the curvature seems to be cancelling the effect of time dilation! 

2.4 Visualizations 

We further understand these results through sketches containing representative 
light-cones of the spacetime. The light-cone structure of the Godel universe as a 
whole in the compact R-cf) plane (see Figure |4]) is much the same as its light-cone 
structure of the infinite r-(j) plane [4J. 

Here too, a CNC forms at the Godel horizon, but once again, does not give 
much insight into how CNGs form in general. This insight is provided by Figure 
El which graphically depicts how the opening and tilting of light-cones allow a 
generic CNC passing through the origin to close. The way this works is this: 
Light-cones beyond the Godel horizon have sections pointing towards the past 
and both dt and d(j) are time-like in this region. Since dt is negative when d0 
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Figure 5: A diagram depicting a Closed Null Curve passing through the origin 
from an interesting perspective. The tangent of the worldline rotates on the 
light-cone, and the light-cones beyond the Godel Horizon have a past-directed 
section as well. 



is positive and vice versa, the signature of the metric is maintained: t moves 
back locally but time does not, because we simultaneously move forward in (jj 
which is also time-like. In what ratio we move back in one and forth in the 
other depends on where in R we are. By going to an appropriately large value 
oi D {D > DcTC, to be precise) the variation in R over the orbit becomes 
sufficiently large that, globally, we end up moving back in both t and (p. This 
is where acausal behaviour sets in. 

Additional insights into the emergence of CNCs/CTCs can be gleaned by 
generalizing the portrayal by Hawking and Ellis of the refocussing of photons 
to WSAs. 

We start by reproducing, in Figure [51 the Hawking and Ellis picture for Rcj) 
space. Since spatial infinity is brought to a finite value, this space is somewhat 
easier to both depict and visualise. 

Since WSAs have all the symmetries of geodesies, WSAs diverging from the 
origin with identical values of v and D (or equivalently. A) will refocus after a 
definite period, T{v, D). In Figure[71 we examine this refocussing for WSAs with 
V = 1 and various values of D (including D = Dj^iuiiGeodisic and D = Dctc) 
by drawing the envelopes of the diverging and reconverging worldlines. 

We note for D = Dctc, the WSAs refocus at the space-time event of their 
departure. In each of the envelopes, we have traced a representative WSA. 
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ClnStUl Null QUPJ* 



Figure 6: The Hawking-EUis picture in {R, (f) space: Envelopes of the null 
geodesies starting from the origin at the same t in different initial directions 0o- 
A cusp develops at the Godel horizon, which itself forms a Closed Null Curve. 
This diagram should be compared to its counterpart, Figure llll in Kundt's 
coordinates. 
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Figure 7: Envelopes of all the v = 1 worldlines emanating from the origin at 

the same t in different initial directions (j)Q, for various fixed vahies of D (or, 
equivalently, of A). The photons refocus at the origin. For D = I/V2, that 
is, for null geodesies, the envelope has a cusp at the Godel horizon. Note that 
the apex of the (lower) cone of divergence approaches that of the (upper) cone 
of reconvergence as D increases. For D = 2-\/2/3, the value corresponding to 
CNCs, the photons reconverge at the same point in both space and time. On 
each of the 5 envelopes, we display a representative (null) worldline. 
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Figure 8: The i(0)-plots of u = 1 worldlines for various values of D. (Note: 
tmid is the values of t at = 0o + f)- To get the curve in {t, i?, 0)-space, curl up 
the 0-axis, aligning the (p^ (f>Q = 0, 1 lines with each other, place the cylinder so 
obtained with 4>—(j)a = along the t-axis at the origin, and rescale the cylinder's 
diameter to D. Note that At — t[Tr) — t(0) is a maximum for D = 1/V2, i.e. 
for null geodesies. Thus acceleration, whatever be its sign, brings the photon 
back earlier to the base-station. This interesting property does not carry over 
to material particles (because D geodesic ~ whereas, DmaxAT — l/^V^, 

and the two are the same only for v = 1). 



In Figure [H we cut open out and flatten out the cylinders on which these 
representatives lie. 

We note that t{(j)) is monotonic only for the D = 0.5 and l/\/2 curves, 
which lie entirely within the Godel horizon. For larger values of D, the slope 
of t{(j)) is in part negative. As explained above, this does not mean that we are 
moving backward in time here. Indeed, we are moving forward in cj) which is 
now likewise time-like. 



3 Kinematics - Rectangular Coordinates 

In the previous section, we solved the geodesic equation only for geodesies that 
pass through the origin because the solutions for geodesies that do not are in 
general much 'dirtier', As a consequence, the accelerated worldlines to which we 
extended the geodesies passed through the origin as well. Since the translational 
symmetry of the Godel universe is not manifest in the cylindrical coordinate sys- 
tem, we could not generalize our results straightforwardly to arbitrary geodesies 
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and WSAs. 

In this section, we accordingly migrate our results to the rectangular coordi- 
nate system {x'^, x, y) with y > 0, introduced by Kundt 12J- This has a manifest 
translational symmetry, which we can use to free our worldlines of the need to 
pass through any fixed point. The line clement is: 

The transformation relations between the rectangular and cylindrical coor- 
dinates are 

2i?sin(/) 

= i+i?2+2i?cos0 

Lux^ = Ljf + 2arctan| ^ tan— )—(?!) (42) 

\1 + R 2J ^ ^ ' 

The xy-projections of the WSAs can be obtained in rectangular coordinates 
by eliminating R and (j> from R — Dsm{(j>~ 0o) using (j40| and (j4T|l . The result 
is 



with 



(x - XcY + iy- VcY = [ycDY (43) 



Dcos<j>o 1 
1 — sm ffln 1 — iJ sm 0n 



We see that the circular projections of diameter D passing through the origin 
in the R-(j) plane are likewise circles in the xy plane centred at (xc^yc), with 
radius ycD. The values of Xc and yc are determined completely by D and 0o- 
All these trajectories pass through the point x = 0, y — w"^ which corresponds 
to the origin i? = in cylindrical coordinates. However, since the radius, ycD 
is (/)o-dependent, geodesies emanating in different directions from (0,a;^^) with 
the same D describe circles of different radii (see Figure [9]). 

To free these trajectories from having to pass through the origin, we simply 
use the most general translation in the plane of rotation of the Godel universe 
(cc, y, x^) —>■ {ax + b, ay, x^) with a > 0, which leaves the metric form-invariant, 
and hence the dynamics unchanged. Under this transformation, the worldline 
continues to be given by (|43p . except that Xc — > axc + b and j/c ciyc- By an 
appropriate choice of a and &, Xc can be made to take any value in (— oo, oo) and 
yc any value in (0, oo). Thus the most general worldline of constant acceleration 
and constant speed can be written parametrically as 

X = Xc + {ycD) cos 9 

y = yc + {ycD) sine (45) 

for any Xc ^ TZ and yc > 0! 

Note further that the transformations between x,y and R,(f) do not involve 
t or x'^. This implies that locally, the transformations between {dt, dfi, dci,) and 
{dxO,dx,dy) are free of boosts. Hence the speed of a particle is identical as 
observed by comoving observers of both coordinate systems. In particular, the 
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Figure 9: The x-y projections of the geodesies of f = 1, and ^ diverging 

from and refocussing at {x,y) = (0,w~^), which corresponds to the origin of 
the cyhndrical coordinates. The Godel Horizon corresponding to this point is 
shown. Note that particles at the same speed leaving {0,lo~^) simultaneously 
in different directions describe circles of different radii, but return at the same 
instant of time. This refocussing is a consequence of the underlying rotational 
symmetry. 
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equation of the geodesic of speed v is ps)) with D = v/^/2. The magnitude of 
acceleration is in any case a 4-scalar, and hence invariant. Thus the acceleration 
of a worldline of speed v with trajectory (j45p for a particular value of D is given 
by exactly the same relation as in cylindrical coordinates: 

2VU! 



A 



1 



1 



DV2 



(46) 



Using the differentials of relations (jin| . (j¥T|) and (jH]), we can transform the 
components of the acceleration in cylindrical corrdinates to those in rectangular 
coordinates. The covariant components of the accaleration at the point 9 on 
the worldline ([iS]) are then 

A 

(0, -cos 61, -sin 61) (47) 



and the contravariant components are 



AP = AujyV2 



cos 9 



Loy 



-, cos 9, sin 9 



(48) 



3.1 The X and s integrals 

To complete the description of the worldlines in the rectangular coordinates, we 
need to relate a;*^ and s to the trajectory parameter 9. The speed as given by 
dSD is: 

2 ^ dx^ + dy^ ^ dx^ + dy^ 

^ 2{ijjydx° + dxY 2{ijjydsy + [dx^ + dy"^) ^ ' 

Using (|45|) . these can be rearranged into differential equations for x^{9) and 
s{9): 

dx° D{sm9 + l/vV2) 



d9 
ds 
d9 



1 



Dsm9 
? D 



vV2 1 + Dsm9 



which can be recast as 
w(x° - a;[J) = 

uj{s - So) 



d9-{2- DV2/v 



dtan ■ 



Dy/2 



dtan ': 



(i? + tan|)2 



[1- D^) + {D + tan 



and readily integrated for \D\ < 1 tco 

2 - DV2/v 



ivix'' - x^j) 



w(s - So) 



3)- 



DV2 



arctan 



J So 



arctan 




(50) 
(51) 



(52) 



(53) 



-'^'^In this case, the arctan has been defined with range [-^, + tt). This solution is valid 
for one orbit ie. for 9 S [So, Oo + 2it). To extend the solution to any number of orbits, add 
niT to the arctan, where n is the number of completed orbits. 
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Figure 10: A map of the future lightcones of the Godel universe in the Kundt 
coordinates. All the light cones contain the x and axes. The only variation 
is along y: As we go further out in y, the cones get stretched along the axis 
by a factor proportional to y^^. 



for \D\ = 1 to 



c^(a;°-a;°) = {6 - 9o) - {2 - V2/v) 



4 



D + tan ■ 



Lli{s - So) 

and for \D\ > 1 to 



'-V2 



1 



D + tan ■ 



uj{x° - X°o) = 
w(s-so) = 



DV2/v 



2 dV2 



tanh 



D + tan ■ 



1 



1 



, , / D + tan f 
tanh"^ 2 
V VD^ - 1 



3.2 Visualizations 



(54) 



(55) 



We now capture the above information in a series of diagrams, as we did for the 
cylindrical system. In Figure fTUl we display Kundt's version of this space-time's 
light-cone structure. 

The translational symmetry along the and x directions is manifest, whereas 
the cones for different values of y are simply scaled versions of one another. The 
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Closed Null Curve 
(Godel Horizon) 



Figure 11: The Hawking-EUis picture in Kundt coordinates: Envelopes of the 
null geodesies starting from the origin at the same t in different initial directions 
(pQ. The Godel horizon forms a Closed Null Curve. The diagram should be 
studied in conjunction with the Figures [H] and [51 The underlying rotational 
symmetry reveals itself in the refocussing of the photons. The apex of the cone 
of divergence and that of reconvergence have the same value of x and y (See 
dotted line). The top view can perhaps be better understood in terms of the 
x-y projections given in Figure [HI 

price we seemingly pay for maintaining manifest homogeneity is that rotation- 
ally symmetry all but disappears from the diagram. A vestige of the underlying 
rotation is nevertheless visible in the definite orientation of all the light-cones, 
which confers on all geodesies a definite sense (see Figure [S]). Moreover, even 
though geodesies diverging from a given point with the same value of v end up 
describing circles of different radii, they continue to refocus0 (see Figure [TT|) . 

Note further that all the future light-cones of this coordinate system contain 
both the a:°- and a;-axes. We can therefore move back in x^ without moving 
back locally in time per se, because whenever we do so, we simultaneously move 
forward in x, which is also time-like. In what ratio we advance in one and 
retreat in the other depends on where in y we happen to be. By adjusting D, 
we can make the variation in y (and thereby the relative tilt in the light-cones) 
sufficient over the orbit to allow a CNC/CTC to close. We capture this closure 
for this system in Figure fT2l 

4 Summary 

To summarize, we have studied worldlines corresponding to motion at constant 
speed and magnitude of acceleration (WSAs). Through these, we can access 

^^The analog of Figure [7] is rather convoluted while that of Figure [8] is very strongly (pg- 
dependent. Since neither is very enlightening, we omit them both. 
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any point in the past or future of a comoving observer at the 'base-station' 
where, because of the rotational symmetry of the space-time, these worldhnes 
refocus periodically. We have extended the Hawking & Ellis depiction of this 
refocussing to the case of photons subject to a constant acceleration orthogonal 
to their instantaneous velocities. By adjusting the magnitude of this acceleration 
we can make the photons refocus at the same point in both space and time. 

We have, additionally, studied CNCs/CTCs more general than the ones 
which develop at/beyond the Godel horizon. In particular, we have shown how 
a critical amount of relative tilt in the light-cones of the space-time allows these 
curves to close, whenever the light-cones in question enclose a spatial axis in ad- 
dition to the usual temporal one. The critical amount of tilt required for closure 
detemines the minimal size of the CNC. We have explicitly illustrated this for 
CNCs of both the cylindrical and the Kundt coordinate systems. Since WSAs 
lie on cylinders parallel to the time-axis in both systems, we have additionally 
traced the emergence of CTCs on more easily visualized 2-dimensional plots, 
by simply cutting open the cylinders, flattening them out and examining the 
curves they carry as functions of the acceleration. 

Thirdly, we have explored WSAs in their full generality by letting the accel- 
eration range from — oo to -l-oo. In particular we have additionally examined 
both open trajectories, and trajectories that curve with the opposite sense to 
that of geodesies. 

Next, we illustrated all these effects in manifestly homogeneous representa- 
tion of the Godel universe in Kundt's coordinates. 
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Lastly, we have derived formulae relating events in Godelian space and time 
to the speed and acceleration of point masses/photons which wc would use (in a 
controlled manner) to travel between them. By carrying out the entire discussion 
purely in terms of parameters with a direct physical interpretation, we have, we 
believe, made intuitive several of the otherwise unintuitive properties of the 
Godel space-time. 
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